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Abstract. We calculate the one-loop contributions to the tau slepton masses in the Minimal 

Supersymmetric Standard Model in the effective potential approach. For the majority of pa- 

^T) ', rameter space under study, those corrections are shown to elevate the value of the lightest stau 

mass. 
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Recently a full one-loop charge breaking effective potential was calculated for the case where, 
other than H^ and H2, the fields tl and tr also acquire non-zero vevs [y, / and r respectively. 
Its minimisation, and impact thereof on CCB bounds, was undertaken in ref. 0. In this article 
, we will not worry about CCB, instead simply note that the full dependence on / and r in the 

^j' effective potential enables us to calculate the one-loop stau masses - they are approximated by 

,-S^ ■ the second derivatives of the effective potential |^], a procedure shown to give very accurate 

results [^ , at least for the Higgs masses. Recently this approach was used to calculate two-loop 
K^ \ corrections to the CP-even Higgs boson masses [^]. Following the conventions and notation of 

^ ■ ref. [y , we recall that we consider the MSSM with Yukawa couplings set to zero for the first and 
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second generations, the superpotential of the model being given by 

W = \tH2QtR + XbHiQbR + KHiLtr + fiH2Hi . (1) 

Supersymmetry is softly broken by adding to the potential explicit mass terms for the gauginos 
and scalar partners, and bilinear and trilinear terms similar in form to those present in the 
superpotential above but multiplied by coefficients B and Ai. When the fields Hf, H2, tl and 
tr have vevs v\l\f2^ t>2/\/2, //\/2 and t I\f2 respectively, the tree-level potential is given by 

V^ = ^[vlif + r')+Pr^] - ^{ArV, + fiV2)lr + Umlvl + mlvl + mlf + 

4 y/2 ^ 

mlr^) - B^,v,v, + '^{vl- vl - l^ + 2r^f + f^{vl-vl + ff . (2) 

The one-loop contributions to the effective potential are given, as usual, by 

^^.-E^*^<t(log§-|) . (3) 

where the Mq, are the tree- level masses of particles of spin Sq,, M is the renormalisation scale 
and Ua = (— 1)^'^"(2sq, + \-)CaQa-, with Ca the number of colour degrees of freedom and Qa 



counting the number of particle-antiparticle states. The presence in the potential of vevs carrying 
electric charge complicates the mass matrices considerably, by causing mixing between neutral 
and charged fields. For example, the tau sneutrinos are mixed with the charged higgs fields, 
originating a three by three mass matrix [O. 



1 The effective potential approach 



The one-loop contributions to the stau masses in the effective potential approach (e.p.a.) will 
be given by the second derivatives of AVi with respect to I and r, so that we have 
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with {fj,fj} 



{/,r}, and we have used the fact that 
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as follows trivially from analysing the {/,t} dependence of the tree-level potential (Q). With 
mass matrices as large as six by six, it would seem impossible to find analytical expressions 
for (^). We can however apply the same trick that allowed us to compute the derivatives 
dM'^/dvi in ref. |0], to wit: the squared masses M^ are the solutions A of the eigenvalue equation 
F = det([M^] — A 1) =0. This equation - a simple polynomial of degree n for an n x n mass 
matrix - implicitly determines the A in terms of the vevs present in the theory. Using (|5|) as well 
as the implicit function theorem we obtain 
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A=M2 



The derivatives of the determinant F are very easy to compute - taking into account that F is 
of the form 



F — 2_^ aiia2jazk ■ ■ ■ f^ijk... 



*,j,/c,. 



we get 
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which is to say, the derivative of the determinant of an n x n matrix becomes the sum of the 
determinants of n matrices, each identical to the initial matrix except for a line replaced with 
derivatives of the original coefficients. Along similar lines we see that the second derivatives of 
F would originate v? determinants. The final expressions obtained in this way for (H) are quite 
involved and we present them in the next section as a function of the sparticles' masses - these 
can be computed either analytically (the only masses that are more involved are the neutralino's; 
see, for instance, ref. [^) or numerically. At this point we must recall that the e.p.a. gives only 
an approximation to the physical masses - in fact, from general principles one sees that [H, H] 
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where the rjj(p^) are the inverse propagators of the stau sector, given by Q 
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p^5ij - Ml + %j{p 



(10) 



The matrix Mf- stands for the coefficients of the terms in the effective potential which are 



quadratic in the stau fields, and is given by 



Mt 



'ml + i [2{ml,)' - {ml)'] cos(2/?) 
fh'^ [A-r — fi cot /?) 



m!^ [At — /i cot (3) 
ml + [{ml,)' - {ml)'] cos(2/3)^ 



(11) 



where fhj- represents the tau fermion mass. The parameters in this mass matrix are all one- 
loop renormalised. The primes in the masses indicate that these are the masses in the e.p.a. 
approximation, that is, (m^)' = 52(^1 +^2)/4i ("t-I)' = (5''^+fl'2)(^i +^2)/^ and fh'^ = \rVi/\/2. 
These are related to the physical, "unprimed" masses by 
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Our input parameters are the experimental values of the masses, so we express M?- in terms of 
unprimed quantities, so that M?- = M?- + liij - the matrix Mf- is identical in form to eq. (Il 
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except the "primed" masses are replaced by the physical ones, and we have 
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(13) 



So, we can rewrite the inverse propagators of eq. ( [lO| ) as 
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where we define Aij{p'^) as 
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Thus, the e.p.a. expressed by eq. (^) consists of neglecting the quantity Ajj(p^) - which is 
equivalent to calculating the stau self energies IIjj at zero external momentum and neglecting 
the effects of the gauge boson and tau self energies present in Ilij. The careful comparison 
between the diagrammatic and the effective potential approaches of ref. |Q] reached similar 
conclusions for the higgs' masses. It was shown there that the e.p.a. produces extremely good 
results, the masses thus calculated differing from the "real" masses by very small amounts, at 
most ~ 3 GeV. It is reasonable to expect, then, that the e.p.a. -calculated stau masses will be 
comparably accurate. We can also expect the one-loop mass corrections to be small - the second 
derivatives of the masses with order to {/, r} will always produce coefficients multiplied by the 
smaller couplings, g' , §2 and A,-. 

^We follow the conventions of ref. ^ and absorb a factor of i in the definition of F; furthermore, the 11 are 
understood as being the real part of the one-loop renormalised self-energies. 



2 One-loop stau masses 

We now list all one-loop contributions to the stau mass matrix as per eq. (Q). The sparticle 
masses are well known in the literature and the factors n^ are given in ref. |||] 0, so we list 
only the non-zero second derivatives of the masses. We adopt the convention X^.^, to denote 
the second derivative of the quantity X with respect to {vi,Vj} = {^r} (and recall that these 
derivatives are evaluated at / = r = 0). The simplest of these contributions are those of the first 
and second generation squarks and sleptons and the gauge bosons, given by 

(16) 
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For the stop, we have 
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Ci = mf + -Xfvi - ^{vi-vt) , (19) 

and identical expressions for the second stop, with the substitution M~ ^^ M? . For the sbottom 
we have 
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"^We take this opportunity to correct a misprint in ref. |l|, though: the factor ria for the first and second 
generation sneutrinos is 4, counting two generations and the existence of both neutrinos and anti-neutrinos. 



Again, identical expressions for the second sbottoni with trivial substitutions. For the charginos 
we find 
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For the neutralinos and tau lepton the expressions are more complex. Prom the formulae of 
ref. |y] it is easy to find ^ 
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where the index i runs from 1 to 6, the two last entries being ztArfi/v^ (the tau mass) and 
the coefficients A^o , . . . F^o are given by 
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and their derivatives, 
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Notice that this is the second derivative of the mass, not its squared. 
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For the charged higgses, 
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Similar expressions hold for Af ^ ^ - notice that if we perform a tree-level minimisation of the 



H. 



potential, this second eigenvalue is actually zero, reflecting the presence of a Goldstone boson. 
With a one-loop minimisation the tree-level charged higgs mass matrix produces two non-zero 
eigenvalues, albeit the second is quite small when compared to the first one, and therefore has a 
small impact in the one-loop potential. This second eigenvalue can even be negative 0] and its 



contribution to the potential neglected following the arguments of [IC]. For the tau sneutrino 
we have 
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The reason for the sharing of coefficients between the second derivatives of MJ^^ and M~ is the 
mixing discussed in ref. ||l|. For the pseudo-scalar Higgs the second derivatives are given by 
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where the functions A are given by 

Au{Ml^) = \l{hn-Ml;)[2Bt?A^- ^?{aH-Ml;) - Al{cH-Ml;) 
ArAMl^) = Xl{JH-Mlj[2Bf,^Ar- t^^as-MlJ -Al{cH-MJj^) 
AUMl^) = XliR [^i'' [an - M|J + Al {ch - M]^^) - 23/,^ A^' 
and the denominator Dfj by 
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With the replacement M^ ^^ M^ , we obtain the expressions for the second pseudo-scalar ^. 
For the Higgs scalars, we obtain 
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where the denominator Dh has the expression 

Dh{MImI) = {Mfj - Ml) [{hff - Ml) {jg - Ml) - il] 

and the functions B are given by 
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*Which, like in the case of the charged Higgses, is a Goldstone boson for a tree-level minimisation. 
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The new coefficients in these expressions are 
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and the derivatives of {chAh, fn^gn} are listed in ref. p. For completeness, they are 
dH,i = 7(4 A^ + 5' - gl) vi dH,T = ^Ar eH,i = -^ Ar 
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At 1 /2 /„„N 

gH,i = — y=Ai gH,T = -^g V2 ■ (37) 

If we perform the substitution Mfj <-> M| in equations (|33|)- ( ^5| ) we obtain the second derivatives 
of the mass of the heaviest scalar Higgs. Finally, for the staus we have 
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Again, the expressions for M? are obtained from these with a simple replacement. Because 
supersymmetry is softly broken the supertrace of the squared masses must be field-independent, 
so we can verify these formulae by checking that Strd'^M'^/dxdy = 0. 



3 Numerical results and discussion 

We now apply our results for the one-loop stau masses to a vast MSSM parameter space. In order 
to try to take into account the effects of the particles' mass thresholds in the renormalisation 
running of the theory's parameters, we follow the procedure outlined in refs. [^ and use as input 
parameters ai = 0.01667, 02 = 0.032, as = 0.1, m^ = 2.95 GeV, m,- = 1.75 GeV and rrit = 167.2 
GeV, at the scale Mz- Accordingly we take the DRED value for v'^ = vf + V2 = (250.75GeV)^, 
and use the supersymmetric two-loop /3-functions to evolve all parameters between Mz and the 
gauge unification scale M(/, defined as the point where the couplings ai and 02 meet. At Mfj 
we input the values of the soft parameters and determine /U and B by minimising the one-loop 
MSSM potential at the scale M, defined as the maximum of Mz and the input scalar and 
gaugino masses. For the soft parameters our strategy was to choose a random value Mq, rac 
and Ag and let the gaugino and scalar masses and A parameters vary randomly within a 30% 
interval of those central values, thus obtaining over 15000 points with input soft masses roughly 
in the interval [10, 1000] GeV and -4 < v4g < 4 TeV. Further, we have taken 2.5 < tan/? < 6.5 
and considered both possible signs for the // parameter. We then impose experimental bounds 
on the sparticles' masses from ref. [^, except, obviously, the bounds on the stau masses, as we 
are interested in checking whether the one-loop contributions change their values considerably. 
As it turns out, for our choice of parameter space, after all other experimental cuts have been 
applied the remaining points (over 10000 of them) correspond to stau masses above the current 
experimental bound (81 GeV) for all but a handful of points. The results of this "scan" of 
the MSSM can be seen in figures (|l|)- (^). In fig. (|l|) we plot the mass difference between the 
one-loop and tree-level masses for the lightest stau, against the maximum M of the input soft 
masses and Mz - M is of the order of the largest masses present in A 14 and as such should 
constitute a good choice for renormalisation scale. Several observations about this plot: for 
M smaller than about 200 GeV there is no substantial difference between the tree-level and 
one-loop results. For larger values, though, there are sizeable differences, usually smaller than 
10 GeV. We observe that for the majority of points in the chosen parameter space the mass 
of the lightest stau increases. As expected, the one-loop contributions are small (typically less 
than 5%) but we find they are not negligible. Figure (||) is the analogous of the previous one, 
but looking at the one-loop/tree-level mass difference for the heaviest stau - we see the one-loop 
contributions tend to decrease the mass of the heavier stau, by as much as ~ 20 GeV. Again, 
these contributions are only a few percent of the total mass (typically less than 8%) but not at 
all insignificant. Now, since the lightest CP-even higgs boson is likely to be discovered before 
the staus, it is useful to look at the relationship between AMf and Mh (we use the full one-loop 
higgs mass). In fig. (^) we plot these two quantities one against the other - despite the fact that 
the largest values of AM^ occur (naturally) for higher values of the input soft masses, they do 
not necessarily correspond to large values of M^- 

In conclusion, we computed the one-loop contributions to the stau masses in the effective 
potential approach and showed they are usually quite small, but can nevertheless be sizeable. 
In particular, the lighter stau mass is shown to increase for the majority of input values we 
considered. If the staus are ever discovered and their mass measured accurately, these mass 
differences could be instrumental in narrowing the parameter space of the MSSM. Caution must 
be exercised in reading these results, though: the e.p.a., we repeat, is an approximation to the 
real mass. For instance, as discussed in ref. Q, the resulting mass has a small renormalisation 
scale dependence. We confirmed that fact by changing the value of M, doubling it or reducing 
it to 100 GeV, but the resulting changes in AMf are indeed very small. We also recall that the 
e.p.a. results for the Higgs sector were off by a few GeV, but again from ref. Q we observe the 
e.p.a tends to underestimate the real masses. As such, our conclusion regarding the increase of 
the lightest stau mass should hold. Of course, the validity of the e.pa. can only be established 



by performing the full diagrammatic calculation, and perhaps the fact the e.p.a. is predicting 
measurable differences for Mf is sufficient reason to undertake it. We also observe that an 
e.p.a. calculation in the stop sector should yield larger one-loop contributions for the simple 
reason it produces second derivatives of the masses proportional to Af . This work is now under 
preparation. 
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Figure 1: The maximum M of Mz and the input soft masses versus the mass difference between 
the one- loop and tree-level calculated lightest stau mass, AMf = M-~ °°^ 
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Figure 2: The maximum M of Mz and the input soft masses versus the mass difference between 
the one-loop and tree-level calculated heaviest stau mass. 
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Figure 3: Lightest CP-even higgs boson mass (one-loop) versus the mass difference between the 
one-loop and tree-level calculated lightest stau mass. 



